home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / sgbcon.z / sgbcon
Encoding:
Text File  |  2002-10-03  |  4.3 KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGBBBBCCCCOOOONNNN((((3333SSSS))))                                                          SSSSGGGGBBBBCCCCOOOONNNN((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGBCON - estimate the reciprocal of the condition number of a real
  10.      general band matrix A, in either the 1-norm or the infinity-norm,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND, WORK,
  14.                         IWORK, INFO )
  15.  
  16.          CHARACTER      NORM
  17.  
  18.          INTEGER        INFO, KL, KU, LDAB, N
  19.  
  20.          REAL           ANORM, RCOND
  21.  
  22.          INTEGER        IPIV( * ), IWORK( * )
  23.  
  24.          REAL           AB( LDAB, * ), WORK( * )
  25.  
  26. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  27.      These routines are part of the SCSL Scientific Library and can be loaded
  28.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  29.      directs the linker to use the multi-processor version of the library.
  30.  
  31.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  32.      4 bytes (32 bits). Another version of SCSL is available in which integers
  33.      are 8 bytes (64 bits).  This version allows the user access to larger
  34.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  35.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  36.      only one of the two versions; 4-byte integer and 8-byte integer library
  37.      calls cannot be mixed.
  38.  
  39. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  40.      SGBCON estimates the reciprocal of the condition number of a real general
  41.      band matrix A, in either the 1-norm or the infinity-norm, using the LU
  42.      factorization computed by SGBTRF.
  43.  
  44.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  45.      condition number is computed as
  46.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  47.  
  48.  
  49. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  50.      NORM    (input) CHARACTER*1
  51.              Specifies whether the 1-norm condition number or the infinity-
  52.              norm condition number is required:
  53.              = '1' or 'O':  1-norm;
  54.              = 'I':         Infinity-norm.
  55.  
  56.      N       (input) INTEGER
  57.              The order of the matrix A.  N >= 0.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSGGGGBBBBCCCCOOOONNNN((((3333SSSS))))                                                          SSSSGGGGBBBBCCCCOOOONNNN((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      KL      (input) INTEGER
  75.              The number of subdiagonals within the band of A.  KL >= 0.
  76.  
  77.      KU      (input) INTEGER
  78.              The number of superdiagonals within the band of A.  KU >= 0.
  79.  
  80.      AB      (input) REAL array, dimension (LDAB,N)
  81.              Details of the LU factorization of the band matrix A, as computed
  82.              by SGBTRF.  U is stored as an upper triangular band matrix with
  83.              KL+KU superdiagonals in rows 1 to KL+KU+1, and the multipliers
  84.              used during the factorization are stored in rows KL+KU+2 to
  85.              2*KL+KU+1.
  86.  
  87.      LDAB    (input) INTEGER
  88.              The leading dimension of the array AB.  LDAB >= 2*KL+KU+1.
  89.  
  90.      IPIV    (input) INTEGER array, dimension (N)
  91.              The pivot indices; for 1 <= i <= N, row i of the matrix was
  92.              interchanged with row IPIV(i).
  93.  
  94.      ANORM   (input) REAL
  95.              If NORM = '1' or 'O', the 1-norm of the original matrix A.  If
  96.              NORM = 'I', the infinity-norm of the original matrix A.
  97.  
  98.      RCOND   (output) REAL
  99.              The reciprocal of the condition number of the matrix A, computed
  100.              as RCOND = 1/(norm(A) * norm(inv(A))).
  101.  
  102.      WORK    (workspace) REAL array, dimension (3*N)
  103.  
  104.      IWORK   (workspace) INTEGER array, dimension (N)
  105.  
  106.      INFO    (output) INTEGER
  107.              = 0:  successful exit
  108.              < 0: if INFO = -i, the i-th argument had an illegal value
  109.  
  110. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  111.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  112.  
  113.      This man page is available only online.
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.